Single-Cell Gene-Regulatory Networks of Advanced Symptomatic Atherosclerosis.

BACKGROUND: While our understanding of the single-cell gene expression patterns underlying the transformation of vascular cell types during the progression of atherosclerosis is rapidly improving, the clinical and pathophysiological relevance of these changes remains poorly understood. METHODS: Single-cell RNA sequencing data generated with SmartSeq2 (≈8000 genes/cell) in nearly 19 000 single cells isolated during atherosclerosis progression in Ldlr-/-Apob100/100 mice with human-like plasma lipoproteins and from humans with asymptomatic and symptomatic carotid plaques was clustered into multiple subtypes. For clinical and pathophysiological context, the advanced-stage and symptomatic subtype clusters were integrated with 135 tissue-specific (atherosclerotic aortic wall, mammary artery, liver, skeletal muscle, and visceral and subcutaneous, fat) gene-regulatory networks (GRNs) inferred from 600 coronary artery disease patients in the STARNET (Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task) study. RESULTS: Advanced stages of atherosclerosis progression and symptomatic carotid plaques were largely characterized by 3 smooth muscle cells (SMCs), and 3 macrophage subtype clusters with extracellular matrix organization/osteogenic (SMC), and M1-type proinflammatory/Trem2-high lipid-associated (macrophage) phenotypes. Integrative analysis of these 6 clusters with STARNET revealed significant enrichments of 3 arterial wall GRNs: GRN33 (macrophage), GRN39 (SMC), and GRN122 (macrophage) with major contributions to coronary artery disease heritability and strong associations with clinical scores of coronary atherosclerosis severity (SYNTAX/Duke scores). The presence and pathophysiological relevance of GRN39 were verified in 5 independent RNAseq data sets obtained from the human coronary and aortic artery, and primary SMCs and by targeting its top-key drivers, FRZB and ALCAM, in cultured human vascular SMCs. CONCLUSIONS: By identifying and integrating the most gene-rich single-cell subclusters of atherosclerosis to date with a coronary artery disease framework of GRNs, GRN39 was identified and independently validated as being critical for the transformation of contractile SMCs into an osteogenic phenotype promoting advanced-stage, symptomatic atherosclerosis.

Coronary Artery Disease Risk Variant Dampens the Expression of CALCRL by Reducing HSF Binding to Shear Stress Responsive Enhancer in Endothelial Cells In Vitro.

BACKGROUND: CALCRL (calcitonin receptor-like) protein is an important mediator of the endothelial fluid shear stress response, which is associated with the genetic risk of coronary artery disease. In this study, we functionally characterized the noncoding regulatory elements carrying coronary artery disease that risks single-nucleotide polymorphisms and studied their role in the regulation of CALCRL expression in endothelial cells. METHODS: To functionally characterize the coronary artery disease single-nucleotide polymorphisms harbored around the gene CALCRL, we applied an integrative approach encompassing statistical, transcriptional (RNA-seq), and epigenetic (ATAC-seq, chromatin immunoprecipitation assay-quantitative polymerase chain reaction, and electromobility shift assay) analyses, alongside luciferase reporter assays, and targeted gene and enhancer perturbations (siRNA and clustered regularly interspaced short palindromic repeats/clustered regularly interspaced short palindromic repeat-associated 9) in human aortic endothelial cells. RESULTS: We demonstrate that the regulatory element harboring rs880890 exhibits high enhancer activity and shows significant allelic bias. The A allele was favored over the G allele, particularly under shear stress conditions, mediated through alterations in the HSF1 (heat shock factor 1) motif and binding. CRISPR deletion of rs880890 enhancer resulted in downregulation of CALCRL expression, whereas HSF1 knockdown resulted in a significant decrease in rs880890-enhancer activity and CALCRL expression. A significant decrease in HSF1 binding to the enhancer region in endothelial cells was observed under disturbed flow compared with unidirectional flow. CALCRL knockdown and variant perturbation experiments indicated the role of CALCRL in mediating eNOS (endothelial NO synthase), APLN (apelin), angiopoietin, prostaglandins, and EDN1 (endothelin-1) signaling pathways leading to a decrease in cell proliferation, tube formation, and NO production. CONCLUSIONS: Overall, our results demonstrate the existence of an endothelial-specific HSF (heat shock factor)-regulated transcriptional enhancer that mediates CALCRL expression. A better understanding of CALCRL gene regulation and the role of single-nucleotide polymorphisms in the modulation of CALCRL expression could provide important steps toward understanding the genetic regulation of shear stress signaling responses.

Secreted Protein Profiling of Human Aortic Smooth Muscle Cells Identifies Vascular Disease Associations.

BACKGROUND: Smooth muscle cells (SMCs), which make up the medial layer of arteries, are key cell types involved in cardiovascular disease, the leading cause of mortality and morbidity worldwide. In response to microenvironment alterations, SMCs dedifferentiate from a contractile to a synthetic phenotype characterized by an increased proliferation, migration, production of ECM (extracellular matrix) components, and decreased expression of SMC-specific contractile markers. These phenotypic changes result in vascular remodeling and contribute to the pathogenesis of cardiovascular disease, including coronary artery disease, stroke, hypertension, and aortic aneurysms. Here, we aim to identify the genetic variants that regulate ECM secretion in SMCs and predict the causal proteins associated with vascular disease-related loci identified in genome-wide association studies. METHODS: Using human aortic SMCs from 123 multiancestry healthy heart transplant donors, we collected the serum-free media in which the cells were cultured for 24 hours and conducted liquid chromatography-tandem mass spectrometry-based proteomic analysis of the conditioned media. RESULTS: We measured the abundance of 270 ECM and related proteins. Next, we performed protein quantitative trait locus mapping and identified 20 loci associated with secreted protein abundance in SMCs. We functionally annotated these loci using a colocalization approach. This approach prioritized the genetic variant rs6739323-A at the 2p22.3 locus, which is associated with lower expression of LTBP1 (latent-transforming growth factor beta-binding protein 1) in SMCs and atherosclerosis-prone areas of the aorta, and increased risk for SMC calcification. We found that LTBP1 expression is abundant in SMCs, and its expression at mRNA and protein levels was reduced in unstable and advanced atherosclerotic plaque lesions. CONCLUSIONS: Our results unravel the SMC proteome signature associated with vascular disorders, which may help identify potential therapeutic targets to accelerate the pathway to translation.

Translatome profiling reveals Itih4 as a novel smooth muscle cell-specific gene in atherosclerosis.

AIMS: Vascular smooth muscle cells (SMCs) and their derivatives are key contributors to the development of atherosclerosis. However, studying changes in SMC gene expression in heterogeneous vascular tissues is challenging due to the technical limitations and high cost associated with current approaches. In this paper, we apply Translating Ribosome Affinity Purification sequencing (TRAP-Seq) to profile SMC-specific gene expression directly from tissue. METHODS AND RESULTS: To facilitate SMC-specific translatome analysis, we generated SMCTRAP mice, a transgenic mouse line expressing EGFP-tagged ribosomal protein L10a (EGFP-L10a) under the control of the SMC-specific αSMA promoter. These mice were further crossed with the atherosclerosis model Ldlr-/-, ApoB100/100 to generate SMCTRAP-AS mice and used to profile atherosclerosis-associated SMCs in thoracic aorta samples of 15-month-old SMCTRAP and SMCTRAP-AS mice. Our analysis of SMCTRAP-AS mice showed that EGFP-L10a expression was localized to SMCs in various tissues, including the aortic wall and plaque. The TRAP fraction demonstrated high enrichment of known SMC-specific genes, confirming the specificity of our approach. We identified several genes, including Cemip, Lum, Mfge8, Spp1, and Serpina3, that are known to be involved in atherosclerosis-induced gene expression. Moreover, we identified several novel genes not previously linked to SMCs in atherosclerosis, such as Anxa4, Cd276, Itih4, Myof, Pcdh11x, Rab31, Serpinb6b, Slc35e4, Slc8a3, and Spink5. Among them, we confirmed the SMC-specific expression of Itih4 in atherosclerotic lesions using immunofluorescence staining of mouse aortic roots and spatial transcriptomics of human carotid arteries. Furthermore, our more detailed analysis of Itih4 showed its link to coronary artery disease (CAD) through the colocalization of GWAS, splice-QTL, and protein-QTL signals. CONCLUSIONS: We generated a SMC-specific TRAP mouse line to study atherosclerosis and identified Itih4 as a novel SMC-expressed gene in atherosclerotic plaques, warranting further investigation of its putative function in extracellular matrix stability and genetic evidence of causality.

Secreted protein profiling of human aortic smooth muscle cells identifies vascular disease associations.

Background: Smooth muscle cells (SMCs), which make up the medial layer of arteries, are key cell types involved in cardiovascular diseases (CVD), the leading cause of mortality and morbidity worldwide. In response to microenvironment alterations, SMCs dedifferentiate from a "contractile" to a "synthetic" phenotype characterized by an increased proliferation, migration, production of extracellular matrix (ECM) components, and decreased expression of SMC-specific contractile markers. These phenotypic changes result in vascular remodeling and contribute to the pathogenesis of CVD, including coronary artery disease (CAD), stroke, hypertension, and aortic aneurysms. Here, we aim to identify the genetic variants that regulate ECM secretion in SMCs and predict the causal proteins associated with vascular disease-related loci identified in genome-wide association studies (GWAS). Methods: Using human aortic SMCs from 123 multi-ancestry healthy heart transplant donors, we collected the serum-free media in which the cells were cultured for 24 hours and conducted Liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based proteomic analysis of the conditioned media. Results: We measured the abundance of 270 ECM and related proteins. Next, we performed protein quantitative trait locus mapping (pQTL) and identified 20 loci associated with secreted protein abundance in SMCs. We functionally annotated these loci using a colocalization approach. This approach prioritized the genetic variant rs6739323-A at the 2p22.3 locus, which is associated with lower expression of LTBP1 in SMCs and atherosclerosis-prone areas of the aorta, and increased risk for SMC calcification. We found that LTBP1 expression is abundant in SMCs, and its expression at mRNA and protein levels was reduced in unstable and advanced atherosclerotic plaque lesions. Conclusions: Our results unravel the SMC proteome signature associated with vascular disorders, which may help identify potential therapeutic targets to accelerate the pathway to translation.

Colchicine promotes atherosclerotic plaque stability independently of inflammation.

Atherosclerosis is a chronic inflammatory disease which is driven in part by the aberrant trans -differentiation of vascular smooth muscle cells (SMCs). No therapeutic drug has been shown to reverse detrimental SMC-derived cell phenotypes into protective phenotypes, a hypothesized enabler of plaque regression and improved patient outcome. Herein, we describe a novel function of colchicine in the beneficial modulation of SMC-derived cell phenotype, independent of its conventional anti-inflammatory effects. Using SMC fate mapping in an advanced atherosclerotic lesion model, colchicine induced plaque regression by converting pathogenic SMC-derived macrophage-like and osteoblast-like cells into protective myofibroblast-like cells which thickened, and thereby stabilized, the fibrous cap. This was dependent on Notch3 signaling in SMC-derived plaque cells. These findings may help explain the success of colchicine in clinical trials relative to other anti-inflammatory drugs. Thus, we demonstrate the potential of regulating SMC phenotype in advanced plaque regression through Notch3 signaling, in addition to the canonical anti-inflammatory actions of drugs to treat atherosclerosis.

Genome-wide census of ATF4 binding sites and functional profiling of trait-associated genetic variants overlapping ATF4 binding motifs.

Activating Transcription Factor 4 (ATF4) is an important regulator of gene expression in stress responses and developmental processes in many cell types. Here, we catalogued ATF4 binding sites in the human genome and identified overlaps with trait-associated genetic variants. We probed these genetic variants for allelic regulatory activity using a massively parallel reporter assay (MPRA) in HepG2 hepatoma cells exposed to tunicamycin to induce endoplasmic reticulum stress and ATF4 upregulation. The results revealed that in the majority of cases, the MPRA allelic activity of these SNPs was in agreement with the nucleotide preference seen in the ATF4 binding motif from ChIP-Seq. Luciferase and electrophoretic mobility shift assays in additional cellular models further confirmed ATF4-dependent regulatory effects for the SNPs rs532446 (GADD45A intronic; linked to hematological parameters), rs7011846 (LPL upstream; myocardial infarction), rs2718215 (diastolic blood pressure), rs281758 (psychiatric disorders) and rs6491544 (educational attainment). CRISPR-Cas9 disruption and/or deletion of the regulatory elements harboring rs532446 and rs7011846 led to the downregulation of GADD45A and LPL, respectively. Thus, these SNPs could represent examples of GWAS genetic variants that affect gene expression by altering ATF4-mediated transcriptional activation.

The angiopoietin receptor Tie2 is atheroprotective in arterial endothelium.

Leukocytes and resident cells in the arterial wall contribute to atherosclerosis, especially at sites of disturbed blood flow. Expression of endothelial Tie1 receptor tyrosine kinase is enhanced at these sites, and attenuation of its expression reduces atherosclerotic burden and decreases inflammation. However, Tie2 tyrosine kinase function in atherosclerosis is unknown. Here we provide genetic evidence from humans and from an atherosclerotic mouse model to show that TIE2 is associated with protection from coronary artery disease. We show that deletion of Tie2, or both Tie2 and Tie1, in the arterial endothelium promotes atherosclerosis by increasing Foxo1 nuclear localization, endothelial adhesion molecule expression and accumulation of immune cells. We also show that Tie2 is expressed in a subset of aortic fibroblasts, and its silencing in these cells increases expression of inflammation-related genes. Our findings indicate that unlike Tie1, the Tie2 receptor functions as the dominant endothelial angiopoietin receptor that protects from atherosclerosis.

Cross-tissue omics analysis discovers ten adipose genes encoding secreted proteins in obesity-related non-alcoholic fatty liver disease.

BACKGROUND: Non-alcoholic fatty liver disease (NAFLD) is a fast-growing, underdiagnosed, epidemic. We hypothesise that obesity-related inflammation compromises adipose tissue functions, preventing efficient fat storage, and thus driving ectopic fat accumulation into the liver. METHODS: To identify adipose-based mechanisms and potential serum biomarker candidates (SBCs) for NAFLD, we utilise dual-tissue RNA-sequencing (RNA-seq) data in adipose tissue and liver, paired with histology-based NAFLD diagnosis, from the same individuals in a cohort of obese individuals. We first scan for genes that are differentially expressed (DE) for NAFLD in obese individuals' subcutaneous adipose tissue but not in their liver; encode proteins secreted to serum; and show preferential adipose expression. Then the identified genes are filtered to key adipose-origin NAFLD genes by best subset analysis, knockdown experiments during human preadipocyte differentiation, recombinant protein treatment experiments in human liver HepG2 cells, and genetic analysis. FINDINGS: We discover a set of genes, including 10 SBCs, that may modulate NAFLD pathogenesis by impacting adipose tissue function. Based on best subset analysis, we further follow-up on two SBCs CCDC80 and SOD3 by knockdown in human preadipocytes and subsequent differentiation experiments, which show that they modulate crucial adipogenesis genes, LPL, SREBPF1, and LEP. We also show that treatment of the liver HepG2 cells with the CCDC80 and SOD3 recombinant proteins impacts genes related to steatosis and lipid processing, including PPARA, NFE2L2, and RNF128. Finally, utilizing the adipose NAFLD DE gene cis-regulatory variants associated with serum triglycerides (TGs) in extensive genome-wide association studies (GWASs), we demonstrate a unidirectional effect of serum TGs on NAFLD with Mendelian Randomization (MR) analysis. We also demonstrate that a single SNP regulating one of the SBC genes, rs2845885, produces a significant MR result by itself. This supports the conclusion that genetically regulated adipose expression of the NAFLD DE genes may contribute to NAFLD through changes in serum TG levels. INTERPRETATION: Our results from the dual-tissue transcriptomics screening improve the understanding of obesity-related NAFLD by providing a targeted set of 10 adipose tissue-active genes as new serum biomarker candidates for the currently grossly underdiagnosed fatty liver disease. FUNDING: The work was supported by NIH grants R01HG010505 and R01DK132775. The Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. The KOBS study (J. P.) was supported by the Finnish Diabetes Research Foundation, Kuopio University Hospital Project grant (EVO/VTR grants 2005-2019), and the Academy of Finland grant (Contract no. 138006). This study was funded by the European Research Council under the European Union's Horizon 2020 research and innovation program (Grant No. 802825 to M. U. K.). K. H. P. was funded by the Academy of Finland (grant numbers 272376, 266286, 314383, and 335443), the Finnish Medical Foundation, Gyllenberg Foundation, Novo Nordisk Foundation (grant numbers NNF10OC1013354, NNF17OC0027232, and NNF20OC0060547), Finnish Diabetes Research Foundation, Finnish Foundation for Cardiovascular Research, University of Helsinki, and Helsinki University Hospital and Government Research Funds. I. S. was funded by the Instrumentarium Science Foundation. Personal grants to U. T. A. were received from the Matti and Vappu Maukonen Foundation, Ella och Georg Ehrnrooths Stiftelse and the Finnish Foundation for Cardiovascular Research.

ErbB signaling is a potential therapeutic target for vascular lesions with fibrous component.

Background: Sporadic venous malformation (VM) and angiomatosis of soft tissue (AST) are benign, congenital vascular anomalies affecting venous vasculature. Depending on the size and location of the lesion, symptoms vary from motility disturbances to pain and disfigurement. Due to the high recurrence of the lesions, more effective therapies are needed. Methods: As targeting stromal cells has been an emerging concept in anti-angiogenic therapies, here, by using VM/AST patient samples, RNA-sequencing, cell culture techniques, and a xenograft mouse model, we investigated the crosstalk of endothelial cells (EC) and fibroblasts and its effect on vascular lesion growth. Results: We report, for the first time, the expression and secretion of transforming growth factor A (TGFA) in ECs or intervascular stromal cells in AST and VM lesions. TGFA induced secretion of vascular endothelial growth factor (VEGF-A) in paracrine fashion, and regulated EC proliferation. Oncogenic PIK3CA variant in p.H1047R, a common somatic mutation found in these lesions, increased TGFA expression, enrichment of hallmark hypoxia, and in a mouse xenograft model, lesion size, and vascularization. Treatment with afatinib, a pan-ErbB tyrosine-kinase inhibitor, decreased vascularization and lesion size in a mouse xenograft model with ECs expressing oncogenic PIK3CA p.H1047R variant and fibroblasts. Conclusions: Based on the data, we suggest that targeting of both intervascular stromal cells and ECs is a potential treatment strategy for vascular lesions having a fibrous component. Funding: Academy of Finland, Ella and Georg Ehnrooth foundation, the ERC grants, Sigrid Jusélius Foundation, Finnish Foundation for Cardiovascular Research, Jane and Aatos Erkko Foundation, GeneCellNano Flagship program, and Department of Musculoskeletal and Plastic Surgery, Helsinki University Hospital.

Dissecting the polygenic basis of atherosclerosis via disease-associated cell state signatures.

Coronary artery disease (CAD) is a pandemic disease where up to half of the risk is explained by genetic factors. Advanced insights into the genetic basis of CAD require deeper understanding of the contributions of different cell types, molecular pathways, and genes to disease heritability. Here, we investigate the biological diversity of atherosclerosis-associated cell states and interrogate their contribution to the genetic risk of CAD by using single-cell and bulk RNA sequencing (RNA-seq) of mouse and human lesions. We identified 12 disease-associated cell states that we characterized further by gene set functional profiling, ligand-receptor prediction, and transcription factor inference. Importantly, Vcam1+ smooth muscle cell state genes contributed most to SNP-based heritability of CAD. In line with this, genetic variants near smooth muscle cell state genes and regulatory elements explained the largest fraction of CAD-risk variance between individuals. Using this information for variant prioritization, we derived a hybrid polygenic risk score (PRS) that demonstrated improved performance over a classical PRS. Our results provide insights into the biological mechanisms associated with CAD risk, which could make a promising contribution to precision medicine and tailored therapeutic interventions in the future.

Deletion of the murine ortholog of human 9p21.3 locus promotes atherosclerosis by increasing macrophage proinflammatory activity.

Background: Several genome-wide association studies have reported a risk locus for coronary artery disease (CAD) in the 9p21. 3 chromosomal region. This region encodes a lncRNA in the INK4 locus (ANRIL) and its genetic variance has a strong association with CAD, but its mechanisms in atherogenesis remain unclear. Objectives: This study aimed to investigate the role of the murine ortholog of human 9p21.3 locus in atherogenesis in hypercholesterolemic mice. Methods: Murine 9p21.3 ortholog knockout mice (Chr4Δ70kb/Δ70kb ) were crossbred with Ldlr -/- ApoB 100/100 mice, and atherosclerotic plaque size and morphology were analyzed on a standard or a high-fat diet (HFD). The hematopoietic cell-specific effect of Chr4Δ70kb/Δ70kb on atherosclerotic plaque development was studied via bone marrow (BM) transplantation, where Chr4Δ70kb/Δ70kb or wild-type BM was transplanted into Ldlr -/- ApoB 100/100 mice. The role of Chr4Δ70kb/Δ70kb in macrophage M1/M2 polarization was studied. In addition, single-cell sequencing data from human and mouse atheroma were analyzed to show the expression profiles of ANRIL and its murine equivalent, Ak148321, in the plaques. Results: Both systemic and hematopoietic Chr4Δ70kb/Δ70kb increased atherosclerosis in Ldlr -/- ApoB 100/100 mice after 12 weeks of HFD. The systemic Chr4Δ70kb/Δ70kb also elevated the number of circulating leukocytes. Chr4Δ70kb/Δ70kb BMDMs showed enhanced M1 polarization in vitro. Single-cell sequencing data from human and mouse atheroma revealed that ANRIL and Ak148321 were mainly expressed in the immune cells. Conclusion: These data demonstrate that both systemic and BM-specific deletion of the murine 9p21.3 risk locus ortholog promotes atherosclerosis and regulates macrophage pro-inflammatory activity, suggesting the inflammation-driven mechanisms of the risk locus on atherogenesis.

Genetic Regulation of SMC Gene Expression and Splicing Predict Causal CAD Genes.

BACKGROUND: Coronary artery disease (CAD) is the leading cause of death worldwide. Recent meta-analyses of genome-wide association studies have identified over 175 loci associated with CAD. The majority of these loci are in noncoding regions and are predicted to regulate gene expression. Given that vascular smooth muscle cells (SMCs) play critical roles in the development and progression of CAD, we aimed to identify the subset of the CAD loci associated with the regulation of transcription in distinct SMC phenotypes. METHODS: We measured gene expression in SMCs isolated from the ascending aortas of 151 heart transplant donors of various genetic ancestries in quiescent or proliferative conditions and calculated the association of their expression and splicing with ~6.3 million imputed single-nucleotide polymorphism markers across the genome. RESULTS: We identified 4910 expression and 4412 splicing quantitative trait loci (sQTLs) representing regions of the genome associated with transcript abundance and splicing. A total of 3660 expression quantitative trait loci (eQTLs) had not been observed in the publicly available Genotype-Tissue Expression dataset. Further, 29 and 880 eQTLs were SMC-specific and sex-biased, respectively. We made these results available for public query on a user-friendly website. To identify the effector transcript(s) regulated by CAD loci, we used 4 distinct colocalization approaches. We identified 84 eQTL and 164 sQTL that colocalized with CAD loci, highlighting the importance of genetic regulation of mRNA splicing as a molecular mechanism for CAD genetic risk. Notably, 20% and 35% of the eQTLs were unique to quiescent or proliferative SMCs, respectively. One CAD locus colocalized with a sex-specific eQTL (TERF2IP), and another locus colocalized with SMC-specific eQTL (ALKBH8). The most significantly associated CAD locus, 9p21, was an sQTL for the long noncoding RNA CDKN2B-AS1, also known as ANRIL, in proliferative SMCs. CONCLUSIONS: Collectively, our results provide evidence for the molecular mechanisms of genetic susceptibility to CAD in distinct SMC phenotypes.

ApoA-I Nanotherapy Rescues Postischemic Vascular Maladaptation by Modulating Endothelial Cell and Macrophage Phenotypes in Type 2 Diabetic Mice.

BACKGROUND: Diabetes is a major risk factor for peripheral arterial disease. Clinical and preclinical studies suggest an impaired collateral remodeling and angiogenesis in response to atherosclerotic arterial occlusion in diabetic conditions, although the underlying mechanisms are poorly understood. OBJECTIVE: To clarify the cellular and molecular mechanisms underlying impaired postischemic adaptive vascular responses and to evaluate rHDL (reconstituted HDL)-ApoA-I nanotherapy to rescue the defect in type 2 diabetic mouse model of hindlimb ischemia. METHODS AND RESULTS: Hindlimb ischemia was induced by unilateral femoral artery ligation. Collateral and capillary parameters together with blood flow recovery were analyzed from normoxic adductor and ischemic gastrocnemius muscles, respectively, at day 3 and 7 post-ligation. In response to femoral artery ligation, collateral lumen area was significantly reduced in normoxic adductor muscles. Distally, ischemic gastrocnemius muscles displayed impaired perfusion recovery and angiogenesis paralleled with persistent inflammation. Muscle-specific mRNA sequencing revealed differential expression of genes critical for smooth muscle proliferation and sprouting angiogenesis in normoxic adductor and ischemic gastrocnemius, respectively, at day 7 post-ligation. Genes typical for macrophage (Mϕ) subsets were differentially expressed across both muscle types. Cell-specific gene expression, flow cytometry, and immunohistochemistry revealed persistent IFN-I response gene upregulation in arterial endothelial cells, ECs and Mϕs from T2DM mice associated with impaired collateral remodeling, angiogenesis and perfusion recovery. Furthermore, rHDL nanotherapy rescued impaired collateral remodeling and angiogenesis through dampening EC and Mϕ inflammation in T2DM mice. CONCLUSIONS: Our results suggest that an impaired collateral remodeling and sprouting angiogenesis in T2DM mice is associated with persistent IFN-I response in ECs and Mϕs. Dampening persistent inflammation and skewing ECs and Mϕ phenotype toward less inflammatory ones using rHDL nanotherapy may serve as a potential therapeutic target for T2DM peripheral arterial disease.

Corrigendum: Exploiting Glutamine Consumption in Atherosclerotic Lesions by Positron Emission Tomography Tracer (2S,4R)-4-18F-Fluoroglutamine.

[This corrects the article DOI: 10.3389/fimmu.2022.821423.].

The Power of Single-Cell RNA Sequencing in eQTL Discovery.

Genome-wide association studies have successfully mapped thousands of loci associated with complex traits. During the last decade, functional genomics approaches combining genotype information with bulk RNA-sequencing data have identified genes regulated by GWAS loci through expression quantitative trait locus (eQTL) analysis. Single-cell RNA-Sequencing (scRNA-Seq) technologies have created new exciting opportunities for spatiotemporal assessment of changes in gene expression at the single-cell level in complex and inherited conditions. A growing number of studies have demonstrated the power of scRNA-Seq in eQTL mapping across different cell types, developmental stages and stimuli that could be obscured when using bulk RNA-Seq methods. In this review, we outline the methodological principles, advantages, limitations and the future experimental and analytical considerations of single-cell eQTL studies. We look forward to the explosion of single-cell eQTL studies applied to large-scale population genetics to take us one step closer to understanding the molecular mechanisms of disease.

Exploiting Glutamine Consumption in Atherosclerotic Lesions by Positron Emission Tomography Tracer (2S,4R)-4-18F-Fluoroglutamine.

Increased glutamine metabolism by macrophages is associated with development of atherosclerotic lesions. Positron emission tomography/computed tomography (PET/CT) with a glutamine analog (2S,4R)-4-18F-fluoroglutamine (18F-FGln) allows quantification of glutamine consumption in vivo. Here, we investigated uptake of 18F-FGln by atherosclerotic lesions in mice and compared the results with those obtained using the glucose analog 2-deoxy-2-18F-fluoro-D-glucose (18F-FDG). Uptake of 18F-FGln and 18F-FDG by healthy control mice (C57BL/6JRj) and atherosclerotic low-density lipoprotein receptor-deficient mice expressing only apolipoprotein B100 (LDLR-/-ApoB100/100) was investigated. The mice were injected intravenously with 18F-FGln or 18F-FDG for in vivo PET/CT imaging. After sacrifice at 70 minutes post-injection, tracer uptake was analyzed by gamma counting of excised tissues and by autoradiography of aorta cryosections, together with histological and immunohistochemical analyses. We found that myocardial uptake of 18F-FGln was low. PET/CT detected lesions in the aortic arch, with a target-to-background ratio (SUVmax, aortic arch/SUVmean, blood) of 1.95 ± 0.42 (mean ± standard deviation). Gamma counting revealed that aortic uptake of 18F-FGln by LDLR-/-ApoB100/100 mice (standardized uptake value [SUV], 0.35 ± 0.06) was significantly higher than that by healthy controls (0.20 ± 0.08, P = 0.03). More detailed analysis by autoradiography revealed that the plaque-to-healthy vessel wall ratio of 18F-FGln (2.90 ± 0.42) was significantly higher than that of 18F-FDG (1.93 ± 0.22, P = 0.004). Immunohistochemical staining confirmed that 18F-FGln uptake in plaques co-localized with glutamine transporter SLC7A7-positive macrophages. Collectively these data show that the 18F-FGln PET tracer detects inflamed atherosclerotic lesions. Thus, exploiting glutamine consumption using 18F-FGln PET may have translational relevance for studying atherosclerotic inflammation.

Single-Cell RNA-Seq Analysis of Olfactory Mucosal Cells of Alzheimer's Disease Patients.

Olfaction is orchestrated by olfactory mucosal cells located in the upper nasal cavity. Olfactory dysfunction manifests early in several neurodegenerative disorders including Alzheimer's disease, however, disease-related alterations to the olfactory mucosal cells remain poorly described. The aim of this study was to evaluate the olfactory mucosa differences between cognitively healthy individuals and Alzheimer's disease patients. We report increased amyloid-beta secretion in Alzheimer's disease olfactory mucosal cells and detail cell-type-specific gene expression patterns, unveiling 240 differentially expressed disease-associated genes compared to the cognitively healthy controls, and five distinct cell populations. Overall, alterations of RNA and protein metabolism, inflammatory processes, and signal transduction were observed in multiple cell populations, suggesting their role in Alzheimer's disease-related olfactory mucosa pathophysiology. Furthermore, the single-cell RNA-sequencing proposed alterations in gene expression of mitochondrially located genes in AD OM cells, which were verified by functional assays, demonstrating altered mitochondrial respiration and a reduction of ATP production. Our results reveal disease-related changes of olfactory mucosal cells in Alzheimer's disease and demonstrate the utility of single-cell RNA sequencing data for investigating molecular and cellular mechanisms associated with the disease.

Functional noncoding SNPs in human endothelial cells fine-map vascular trait associations.

Functional consequences of genetic variation in the noncoding human genome are difficult to ascertain despite demonstrated associations to common, complex disease traits. To elucidate properties of functional noncoding SNPs with effects in human endothelial cells (ECs), we utilized our previous molecular quantitative trait locus (molQTL) analysis for transcription factor binding, chromatin accessibility, and H3K27 acetylation to nominate a set of likely functional noncoding SNPs. Together with information from genome-wide association studies (GWASs) for vascular disease traits, we tested the ability of 34,344 variants to perturb enhancer function in ECs using the highly multiplexed STARR-seq assay. Of these, 5711 variants validated, whose enriched attributes included: (1) mutations to TF binding motifs for ETS or AP-1 that are regulators of the EC state; (2) location in accessible and H3K27ac-marked EC chromatin; and (3) molQTL associations whereby alleles associate with differences in chromatin accessibility and TF binding across genetically diverse ECs. Next, using pro-inflammatory IL1B as an activator of cell state, we observed robust evidence (>50%) of context-specific SNP effects, underscoring the prevalence of noncoding gene-by-environment (GxE) effects. Lastly, using these cumulative data, we fine-mapped vascular disease loci and highlighted evidence suggesting mechanisms by which noncoding SNPs at two loci affect risk for pulse pressure/large artery stroke and abdominal aortic aneurysm through respective effects on transcriptional regulation of POU4F1 and LDAH Together, we highlight the attributes and context dependence of functional noncoding SNPs and provide new mechanisms underlying vascular disease risk.

Meta-Analysis of Smooth Muscle Lineage Transcriptomes in Atherosclerosis and Their Relationships to In Vitro Models.

BACKGROUND: Vascular smooth muscle cells (VSMC) exhibit phenotypic plasticity in atherosclerotic plaques, and among other approaches, has been modeled in vitro by cholesterol loading. METHODS: Meta-analysis of scRNA-seq data from VSMC lineage traced cells across five experiments of murine atherosclerosis was performed. In vivo expression profiles were compared to three in vitro datasets of VSMCs loaded with cholesterol and three datasets of polarized macrophages. RESULTS: We identified 24 cell clusters in the meta-analysis of single cells from mouse atherosclerotic lesions with notable heterogeneity across studies, especially for macrophage populations. Trajectory analysis of VSMC lineage positive cells revealed several possible paths of state transitions with one traversing from contractile VSMC to macrophages by way of a proliferative cell cluster. Transcriptome comparisons between in vivo and in vitro states underscored that data from three in vitro cholesterol-treated VSMC experiments did not mirror cell state transitions observed in vivo. However, all in vitro macrophage profiles analyzed (M1, M2, and oxLDL) were more similar to in vivo profiles of macrophages than in vitro VSMCs were to in vivo profiles of VSMCs. oxLDL loaded macrophages showed the most similarity to in vivo states. In contrast to the in vitro data, comparison between mouse and human in vivo data showed many similarities. CONCLUSIONS: Identification of the sources of variation across single cell datasets in atherosclerosis will be an important step towards understanding VSMC fate transitions in vivo. Also, we conclude that cholesterol-loading in vitro is insufficient to model the VSMC cell state transitions observed in vivo, which underscores the need to develop better cell models. Mouse models, however, appear to reproduce a number of the features of VSMCs in human plaques.